Genome-wide classification of epigenetic activity reveals regions of enriched heritability in immune-related traits.

Epigenetics underpins the regulation of genes known to play a key role in the adaptive and innate immune system (AIIS). We developed a method, EpiNN, that leverages epigenetic data to detect AIIS-relevant genomic regions and used it to detect 2,765 putative AIIS loci. Experimental validation of one of these loci, DNMT1, provided evidence for a novel AIIS-specific transcription start site. We built a genome-wide AIIS annotation and used linkage disequilibrium (LD) score regression to test whether it predicts regional heritability using association statistics for 176 traits. We detected significant heritability effects (average |τ∗|=1.65) for 20 out of 26 immune-relevant traits. In a meta-analysis, immune-relevant traits and diseases were 4.45× more enriched for heritability than other traits. The EpiNN annotation was also depleted of trans-ancestry genetic correlation, indicating ancestry-specific effects. These results underscore the effectiveness of leveraging supervised learning algorithms and epigenetic data to detect loci implicated in specific classes of traits and diseases.

Prior fluid and electrolyte imbalance is associated with COVID-19 mortality.

Background: The COVID-19 pandemic represents a major public health threat. Risk of death from the infection is associated with age and pre-existing comorbidities such as diabetes, dementia, cancer, and impairment of immunological, hepatic or renal function. It remains incompletely understood why some patients survive the disease, while others do not. As such, we sought to identify novel prognostic factors for COVID-19 mortality. Methods: We performed an unbiased, observational retrospective analysis of real world data. Our multivariable and univariable analyses make use of U.S. electronic health records from 122,250 COVID-19 patients in the early stages of the pandemic. Results: Here we show that a priori diagnoses of fluid, pH and electrolyte imbalance during the year preceding the infection are associated with an increased risk of death independently of age and prior renal comorbidities. Conclusions: We propose that future interventional studies should investigate whether the risk of death can be alleviated by diligent and personalized management of the fluid and electrolyte balance of at-risk individuals during and before COVID-19.

A mega-analysis of expression quantitative trait loci (eQTL) provides insight into the regulatory architecture of gene expression variation in liver.

Genome-wide association studies (GWAS) have identified numerous genetic variants in the human genome associated with diseases and traits. Nevertheless, for most loci the causative variant is still unknown. Expression quantitative trait loci (eQTL) in disease relevant tissues is an excellent approach to correlate genetic association with gene expression. While liver is the primary site of gene transcription for two pathways relevant to age-related macular degeneration (AMD), namely the complement system and cholesterol metabolism, we explored the contribution of AMD associated variants to modulate liver gene expression. We extracted publicly available data and computed the largest eQTL data set for liver tissue to date. Genotypes and expression data from all studies underwent rigorous quality control. Subsequently, Matrix eQTL was used to identify significant local eQTL. In total, liver samples from 588 individuals revealed 202,489 significant eQTL variants affecting 1,959 genes (Q-Value < 0.001). In addition, a further 101 independent eQTL signals were identified in 93 of the 1,959 eQTL genes. Importantly, our results independently reinforce the notion that high density lipoprotein metabolism plays a role in AMD pathogenesis. Taken together, our study generated a first comprehensive map reflecting the genetic regulatory landscape of gene expression in liver.

The role of noncoding RNAs in chromatin regulation during differentiation.

A myriad of nuclear noncoding RNAs (ncRNAs) have been discovered since the paradigm of RNAs as plain conveyors of protein translation was discarded. There is increasing evidence that at vital intersections of developmental pathways, ncRNAs target the chromatin modulating machinery to its site of action. However, the mechanistic details of processes involved are still largely unclear, and well-characterized metazoan ncRNA species implicated in chromatin regulation during differentiation remain few. Nevertheless, four major categories are slowly emerging: cis-acting antisense ncRNAs that flag the neighboring genes for the propagation of chromatin marks; allele-specific ncRNAs that perform similar tasks, but target larger loci that typically vary in size from hundreds of thousands of base pairs to a whole chromosome; structural ncRNAs proposed to act as scaffolds that couple chromatin shaping complexes of distinct functionalities; and cofactor ncRNAs with a capacity to inhibit or activate essential components of the intertwined chromatin and transcription apparatuses.

Nuclear-localized tiny RNAs are associated with transcription initiation and splice sites in metazoans.

We have recently shown that transcription initiation RNAs (tiRNAs) are derived from sequences immediately downstream of transcription start sites. Here, using cytoplasmic and nuclear small RNA high-throughput sequencing datasets, we report the identification of a second class of nuclear-specific approximately 17- to 18-nucleotide small RNAs whose 3' ends map precisely to the splice donor site of internal exons in animals. These splice-site RNAs (spliRNAs) are associated with highly expressed genes and show evidence of developmental stage- and region-specific expression. We also show that tiRNAs are localized to the nucleus, are enriched at chromatin marks associated with transcription initiation and possess a 3'-nucleotide bias. Additionally, we find that microRNA-offset RNAs (moRNAs), the miR-15/16 cluster previously linked to oncosuppression and most small nucleolar RNA (snoRNA)-derived small RNAs (sdRNAs) are enriched in the nucleus, whereas most miRNAs and two H/ACA sdRNAs are cytoplasmically enriched. We propose that nuclear-localized tiny RNAs are involved in the epigenetic regulation of gene expression.

Nucleosomes are preferentially positioned at exons in somatic and sperm cells.

Nucleosome positioning is constrained at eukaryotic transcription start sites and implicated in transcriptional regulation. Moreover, recent observations indicate that chromatin structure, transcription and splicing are functionally intertwined, and that modified nucleosomes with trimethylation of lysine 36 in histone subunit 3 (H3K36me3) are enriched at internal exons and the downstream flanking intronic regions of highly expressed genes. However, the position of nucleosomes in the interior of genes has been thought to be largely random. Here we show, by analysis of data sets from human sperm and T cells and medaka (Japanese killifish, Oryzias latipes) blastulae, that internal exons of genes are characterized by sharply elevated average nucleosome occupancy in comparison to flanking intronic sequences. We also show that the preferential positioning of nucleosomes at internal exons is independent of their modification status, and of the GC content, conservation or the expression level of the exon. These findings show that the location of exons is recorded in the chromatin structure and may be inherited across generations. Such embedded information may underpin transcriptionally coupled exon recognition and splice site selection.

Molecular evolution of the HBII-52 snoRNA cluster.

HBII-52 is a human brain-specific C/D box snoRNA that potentially regulates the editing and/or alternative splicing of the serotonin receptor. Forty-two nearly identical copies of the HBII-52 gene are located immediately downstream of the SNRPN protein-coding gene in an imprinted locus associated with Prader-Willi syndrome. Other eutherian mammals, with genomic assemblies covering the corresponding locus, also have multiple orthologous copies of HBII-52. The SNRPB gene, which is known to have given rise to SNRPN through gene duplication, expresses a C/D box snoRNA, SNORD119, from its fifth intron. Here we show that, despite the fact that they lie in different positions relative to the orthologous SNRPB/SNRPN coding sequences, there are significant sequence similarities between SNORD119 and HBII-52, including the antisense element and the stem-forming regions. By analysing these snoRNAs in marsupial and eutherian mammal genomes, we reconstruct the likely evolutionary history of the HBII-52 cluster and SNORD119 and suggest that they have evolved from a common ancestor.

Evidence for control of splicing by alternative RNA secondary structures in Dipteran homothorax pre-mRNA.

In a recent study that identified highly evolutionary conserved sequences in three genomes of Diptera species we described an ultraconserved element found at an internal exon-intron junction of the Drosophila melanogaster homothorax (hth) gene that appeared to be involved in the control of hth pre-mRNA splicing. We also discussed a possible role of RNA secondary structure at this site in the regulation of hth pre-mRNA splicing. In this report we identify a shorter evolutionary conserved intronic element within the hth gene that is located downstream of the first element and has sequence complementarity to it. We demonstrate that intramolecular interactions between these two elements would give rise to alternative RNA secondary structures, which in turn may result in differential control of homothorax pre-mRNA splicing. We also provide additional comparative genomic data from several newly available insect genomes supporting our original conclusion that these conserved elements are important in the post-transcriptional regulation of homothorax gene expression in Diptera.